An allosteric model of circadian KaiC phosphorylation 
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In a recent series of ground-breaking experiments, Nakajima et al. [Science 308, 414-415 (2005)] 
showed that the three cyanobacterial clock proteins KaiA, KaiB, and KaiC are sufficient in vitro 
to generate circadian phosphorylation of KaiC. Here, we present a mathematical model of the Kai 
system. At its heart is the assumption that KaiC can exist in two conformational states, one 
favoring phosphorylation and the other dephosphorylation. Each individual KaiC hexamer then has 
a propensity to be phosphorylated in a cyclic manner. To generate macroscopic oscillations, however, 
the phosphorylation cycles of the different hexamers must be synchronized. We propose a novel 
synchronisation mechanism based on differential affinity: KaiA stimulates KaiC phosphorylation, 
but the limited supply of KaiA dimers binds preferentially to those KaiC hexamers that are falling 
behind in the oscillation. KaiB sequesters KaiA and stabilizes the dephosphorylating KaiC state. 
We show that our model can reproduce a wide range of published data, including the observed 
insensitivity of the oscillation period to variations in temperature, and that it makes nontrivial 
predictions about the effects of varying the concentrations of the Kai proteins. 



Cyanobacteria are the simplest organisms to use circa- 
dian rhythms to anticipate the changes between day and 
night. In the cyanobacterium Synechococcus elongatus, 
the three genes kai A, kaiB and fcazCare the central com- 
ponents of the circadian clock [l| . In higher organisms, it 
is believed that circadian rhythms arc driven primarily by 
transcriptional feedback^. KaiC phosphorylation, how- 
ever, shows a circadian rhythm even when transcription 
and translation are inhibited 0]. Still more remarkably, 
it was recently shown that this rhythmic KaiC phospho- 
rylation can be reconstituted in vitro in the presence of 
only KaiA, KaiB, and ATP [3|. The Kai system thus 
represents a very rare example of a functional biochem- 
ical circuit that can be re-created in the test tube. It is 
a major open question to explain how stable oscillations 
can result from the experimentally observed interactions 
among the different Kai proteins. 

In living cells, KaiC phosphorylation increases during 
the subjective day and decreases during the subjective 
night, and this phosphorylation in turn regulates KaiC's 
activity as a global transcriptional repressor Ml. KaiC 
forms a hexamer both in vivo and in vitro [8|; KaiA 
is present in the cell as a dimer Q and KaiB as a 
dimer [H, Q or a tetramer [ll|. KaiC has both auto- 
dephosphorylation and weaker auto-phosphorylation ac- 
tivity, with the latter dependent on ATP binding 0, || 
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KaiC phosphorylation is stimulated by 
whereas KaiB appears to interfere with 
121 . Ii~3l . KaiC hexamers form hetero- 
multimeric complexes with KaiA and KaiB dimers, but 
one such complex contains no more than one KaiC hex- 
amer [H, H, El]. The composition of these complexes 
varies with a roughly 24 hour period. 

Nakajima et aL's striking observation of in vitro os- 
cillations [H in KaiC phosphorylation poses an obvious 
challenge for modelers. Not only is there the potential for 
detailed comparisons between a model's predictions and 
the wealth of experimental data, the Kai system also has 
several novel features. Most notably, ATP is consumed, 
and the system is driven out of equilibrium, only through 
the repeated phosphorylation and dephosphorylation of 
KaiC. Other reactions, such as the (un)binding of KaiA 
and KaiB to KaiC, should thus obey detailed balance. 
Moreover, unlike in most biological oscillations in 
the Kai system the proteins are neither created nor de- 
stroyed. This imposes significant constraints on any 
model that hopes to explain the in vitro oscillations. 

Several previous papers have put forward interesting 
ideas on how these oscillations might occur [2(| H3] • 
However, they either require that KaiC hexamers can 
bind to each other to form higher order complexes |2fj, 
l2lj ]. a possibility ruled out by recent experiments 0, Q 
, or they assume that KaiA and KaiB can each take on 
multiple forms [22]] . In the latter case, the authors pro- 
pose that these forms may correspond to different sub- 
cellular localizations, but that suggestion cannot hold for 
the in vitro system. Embcrly and Wingreen [2p| intro- 
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duced the elegant hypothesis that exchange of monomers 
among KaiC hexamers might contribute to oscillations, 
an idea supported by recent observations [§]. Their own 
work, however, shows that such exchange by itself is in- 
sufficient to produce sustained oscillations. There is thus 
clearly another mechanism at work in the Kai system. 

Here, we propose such a mechanism. Our model is 
built upon two key elements. First, we hypothesize that 
an isolated KaiC hexamer already has a tendency to 
be cyclically phosphorylated and dcphosphorylated as it 
flips between two allosteric states. Second, we suggest 
that these noisy oscillations of individual hexamers can 
be synchronized through the phenomenon of differential 
affinity, whereby the laggards in a population outcom- 
pete the other hexamers for a limited number of KaiA 
molecules that stimulate phosphorylation. The slowest 
hexamers thus speed up while the fastest are forced to 
slow down, causing the entire population to oscillate in 
phase. 

In the rest of this paper, we first, in section I, show how 
a simple picture of allosteric transitions in KaiC leads 
each hexamer to have an intrinsic phosphorylation cycle. 
We then use an idealized model to introduce the concept 
of differential affinity in section II. This model shows that 
the mechanism requires only a few generic ingredients, 
suggesting that the same synchronization principle could 
be at work in other biological systems. Finally, we turn 
in section III to a more complicated model of the Kai 
system. This model reproduces the phosphorylation be- 
havior of KaiC not only in the in vitro experiments in 
which all three Kai proteins are present, but also in sys- 
tems where KaiA and/or KaiB are absent. In fact, we 
found that the experiments on the various subsets of the 
three Kai proteins strongly constrain the model's design. 
Beyond synchronizing oscillations, KaiA and KaiB must 
also bind to and stabilize one or the other KaiC state. 
When this binding is strong enough, the system more- 
over exhibits temperature compensation, as observed 
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FIG. 1: (A) Schematic free energy levels for KaiC subunits. 
Subunits can be in the active (A) or the inactive (I) state. 
Furthermore, subunits can be phosphorylated (p) and bind 
ATP or ADP. Phosphorylation favors the inactive state, nu- 
cleotide binding the active state. (B) Reaction network for a 
KaiC hexamer with 6 phosphorylation sites. C; and d denote 
a hexamer with i phosphorylated monomers in respectively 
the active and inactive state. (C) Phosphorylation cycles for 
the model in B. The phosphorylation level p of a single KaiC 
hexamer, as obtained by stochastic simulations (solid line), 
and of a population of hexamers, obtained from the mean- 
field rate equations (dashed line). The phosphorylation level 

P = E, *([C 4 ] + [Ci])/ E, 6([0 4 ] + [Ci]), where [C»] is the con- 
centration of hexamers in state Ci. For parameter values, see 
Supporting Information. 



I. ALLOSTERIC MODEL 

In this section, we introduce a simple model of al- 
losteric transitions in KaiC that naturally gives rise to re- 
peated rounds of phosphorylation and dephosphorylation 
within each hexamer. Allosteric conformational changes 
are widespread in biochemistry, and the conformations of 
members of the RecA/DnaB superfamily, to which KaiC 
belongs, have been extensively studied [23l. |24|. 

KaiC monomers Although there is strong evidence 
that KaiC monomers can be phosphorylated at multi- 
ple sites 0, HH, [26| , most published data do not distin- 
guish between different phosphorylated forms. We thus 
assume that KaiC monomers have only two phosphory- 
lation states, phosphorylated and unphosphorylated. 

We postulate that an individual KaiC monomer can 
be in either an active (A) or an inactive (I) conforma- 
tion. Fig. [T]A_ shows the free energies of the different 



monomer states; we consider ATP binding only to un- 
phosphorylated and ADP binding only to phosphorylated 
monomers. As the figure indicates, we assume that phos- 
phorylation favors the inactive over the active state. Nu- 
cleotides have a higher affinity for monomers in the active 
state than for those in the inactive state, so nucleotide 
binding favors the active state over the inactive one. We 
also take both the transfer of a phosphate from ATP 
to a KaiC monomer and the removal of the phosphate 
from the monomer to be thcrmodynamically favorable. 
Taken together, these elements allow for a phosphoryla- 
tion cycle: unphosphorylated monomers prefer to be in 
the A state, where ATP hydrolysis drives phosphoryla- 
tion, while phosphorylated monomers prefer to be in the 
I state, where dephosphorylation occurs spontaneously. 
Each monomer thus tends to go through the sequence of 
reactions A -> A- ATP -> Ap-ADP -> Ip-ADP -> Ip -> 
I — ► A, during which one ATP molecule is hydrolyzed. 
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KaiC hexamers. In the spirit of the MWC 
model [27] . we assume that the energetic cost of hav- 
ing two different monomer conformations in the same 
hexamer is prohibitively large. We can then speak of a 
hexamer as being in either the A or the I state. The 
total (free) energy of the hexamer is simply the sum of 
the contributions from its constituent monomers. Highly 
phosphorylatcd hexamers thus prefer to be in the I state, 
where they will be dephosphorylated, while weakly phos- 
phorylated hexamers prefer the A state, where they will 
be phosporylated. As a result, each hexamer tends to go 
through a cycle in which it is first phosphorylated, then 
dephosphorylated, as indicated in Fig. Q}3 and Fig. 2 of 
the Supporting Information. 

The transition (or flip) rates /, for a hexamer with 
i phosphorylatcd monomers to go from the A to the I 
state and bi to go from the I to the A state depend on 
the energy barriers to the conformational changes. If we 
assume that ATP and ADP exchange are fast, so that 
the free energy of each state is well-defined, then the 
difference in free energy AG between the I and A states 
grows linearly with i: AG(z) = iAG p + (6 — i)AG u , where 
the subscripts p and u refer to the free-energy differences 
for phosphorylated and unphosphorylated monomers, re- 
spectively. The natural phenomenological assumption is 
then that the flip rates depend exponentially on the free- 
energy difference: 

fi = k e X p[AG(i)/2) ~ c * (1) 
h = fc exp[-AG(i)/2]^c- 1 , (2) 

where fco sets the basic timescale and c = exp[(AG p — 
AG u )/2]. Alternatively, one can develop an explicit tran- 
sition state theory that includes the number of bound 
nucleotides as one of the order parameters for the confor- 
mational transition (see Supporting Information). This 
leads to flipping rates that vary exponentially with i just 
as in Eqs.HHU In either case, the rates depend strongly 
on the phosphorylation level, with the consequence that 
hexamers can flip from A to I only when most of their 
monomers are phosphorylatcd and from I to A only when 
most are not phosphorylated. 

In Fig. [Tp, we show the time dependence of the 
phosphorylation level of a single KaiC hexamer obtained 
by Monte Carlo simulations of the chemical master 
equation (see Supporting Information) [28| . Initially, the 
KaiC hexamer is in the unphosphorylated active state, 
Go- KaiC (de)phosphorylation clearly occurs in a cyclic 
fashion, with few transitions from one conformation to 
the other occurring at intermediate phosphorylation. 
However, both the amplitude and the period of the 
phosphorylation cycle are highly variable. Due to this 
variability, the phosphorylation cycles of a population of 
independent KaiC hexamers will quickly dephase. As a 
result, in Fig. [Tp the mean phosphorylation level of the 
KaiC population calculated by integrating deterministic 
rate equations based on the law of mass action shows no 
oscillatory behavior. In order to explain the oscillations 
observed in the in vitro Kai system, the uncoupled 



phosphorylation cycles of the individual KaiC hexamers 
need to be synchronized. 



II. SYNCHRONISATION WITH 
DIFFERENTIAL AFFINITY 

The natural candidates to link the phosporylation 
states of different KaiC hexamers are the other two Kai 
proteins. Here, we present a simple model in which KaiA 
plays this role by catalyzing phosphorylation in the ac- 
tive state, while KaiB is completely absent. This model 
will allow us to introduce several important ideas without 
the distractions that a more faithful description would 
entail. It shows synchronized limit-cycle oscillations in 
KaiC phosphorylation, provided that the concentration 
of KaiA is sufficiently small and that KaiA binds to 
KaiC with differential affinity: KaiA should bind most 
strongly to weakly phosphorylated KaiC hexamers. Al- 
though here we limit our discussion to a particular model 
inspired by the Kai system, the differential affinity mech- 
anism is also amenable to a more general, abstract for- 
mulation, described in the Supporting Information, which 
also shows that the oscillations arise through a supercrit- 
ical Hopf bifurcation. 

We assume that only a single dimer of KaiA can bind 
to a KaiC hexamer, and we force every hexamer to pro- 
ceed through the states C0-C6 and Cg-Co in order (thus 
neglecting intermediate flips). This yields 

Cfe ?i i>q n ^o\ 

6 - *■ — * ^0 (3) 

C, ^ a C,_! (4) 
d + A «=* AC, ^ C i+1 + A (i ? 6) . (5) 

We use deterministic, mass-action kinetics to model the 
effects of these reactions. Here and Q denote i-fold 
phosphorylatcd KaiC hexamers in the active and inactive 
states, and A denotes a KaiA dimer. Eqs. describe 
the same processes within a single hexamer as the dia- 
gram in Fig.[Tj3, with the exception that phosphorylation 
of the active state now requires KaiA, which associates 
with active KaiC with on and off rates k M and k^ h and 
stimulates phosphorylation with a rate k p { (Eq. [5]) . We 
implement differential affinity by setting k^ h = k^a 1 , 
with a > 1 (see Supporting Information). 

Fig.[5]A. shows the mean phosphorylation level of a pop- 
ulation of KaiC hexamers as a function of time. In con- 
trast to the behavior seen in Fig.[Tp, there are clear oscil- 
lations: The KaiA dimers effectively couple the phospho- 
rylation cycles of the individual KaiC hexamers. During 
the phosphorylation phase of the oscillations, most hex- 
amers are in the active form. In this state they can bind 
KaiA, which stimulates their phosphorylation. The con- 
centration of KaiA, however, is limited; indeed, in this 
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number of subunits phosphorylated 

FIG. 2: Limit cycle oscillations in KaiC phosphorylation for 
the simplified model defined by Eqs. |3][5] (A). Mean phospho- 
rylation level p and normalized concentration of free KaiA 
[A]/[A]t- During the phosphorylation phase [A] drops al- 
most to zero. (B), (C) and (D) KaiA binding at 3 stages of 
the phosphorylation phase, marked by circles in (A). KaiA 
favors the less phosphorylated KaiC hexamers. [C]t, total 
KaiC concentration; ACi, complex of KaiA and Ci. We take 
[A]t/[C]t = 0.02 and initially set [Co] = [C]t; see Supporting 
Information for other parameters. 



part of the cycle the concentration of free KaiA is close 
to zero (Fig. [2]A). This means that the KaiC hexamers 
compete with one another for KaiA. In this competition, 
the complexes with a lower degree of phosphorylation 
have the advantage, because they have a higher affinity 
for KaiA. Hence, during the phosphorylation phase, KaiA 
will be mostly bound to the lagging hexamers. This is 
shown in Fig. 03, where the concentrations [Ci] and [ACi] 
are plotted versus i for three different time points. The 
distributions do not overlap: KaiA has a clear preference 
for the less phosphorylated KaiC hexamers. Since the 
phosphorylation rate depends on the amount of bound 
KaiA, laggards with a low degree of phosphorylation will 
be phosphorylated at a high rate, whereas front runners 
with a high degree of phosphorylation will be unable to 
increase their phosphorylation level further. This is the 
essence of the differential affinity synchronisation mech- 
anism. 



III. FULL MODEL OF THE KAI SYSTEM 

The simple model of the previous section showed how 
differential affinity can synchronize the oscillations of the 
different KaiC hexamers. This model, however, neglects 
KaiB completely and is not consistent with the large 
body of experimental data on the Kai system. Here, we 
present a more refined allosteric model. 



A. The model 



The key ingredients of our model are: 

1. KaiA can bind to the active form of KaiC, stimulating 
KaiC phosphorylation. 

Recent experiments suggest that, in the absence of KaiB, 
KaiA binds as a single dimer to the CII domain of the 
KaiC hexamer [2t|. Since this is the domain containing 
KaiC's phosphorylation site, it seems reasonable that the 
affinity of KaiA might depend on the phosphorylation 
state of KaiC. We thus assume, as before, that a single 
KaiA can bind to the active state of KaiC and that the 
affinity of KaiA for active KaiC decreases as the phos- 
phorylation level increases. 

2. The active state of KaiC is more stable than the in- 
active one. 

The experiments of [1, Q show that in the presence of 
only KaiA, KaiC becomes very highly phosphorylated. 
In the absence of KaiB, KaiC should thus have no ten- 
dency to cyclically phosphorylate and dcphosphorylatc. 
This requires that the active state of KaiC has a lower 
free energy than the inactive one (thus shifting the en- 
ergy levels in Fig. [TJ'V from their symmetric values). 

3. KaiB can bind to the inactive form of KaiC. The 
resulting KaiB-KaiC complex can then bind to and se- 
quester KaiA. 

The phosphorylation behavior of KaiC in the presence 
of KaiB, but not KaiA, is essentially identical to that of 
KaiC in the absence of both KaiA and KaiB [T^.[T3j. This 
observation strongly suggests that KaiB does not directly 
affect phosphorylation and dcphosphorylation rates. We 
propose instead the following functions for KaiB: 1) KaiB 
can increase the stability of the inactive state of KaiC 
by binding to it. This restores the capacity of individ- 
ual KaiC hexamers to sustain phosphorylation cycles. 2) 
Strong binding of KaiA by KaiB associated with the in- 
active KaiC hexamers reduces the concentration of free 
KaiA dimers. This leads to a variant of the differential 
affinity mechanism, which is necessary for synchroniz- 
ing the oscillations of the different KaiC hexamers, as 
we clarify below. Based upon the measured size of the 
heteromultimeric complexes 0, H|, we assume that the 
inactive form of KaiC can bind two KaiB dimers, and 
that B2C4, B2C3, B2C2, and B2C1 can each bind two 
KaiA dimers with high affinity. Neither assumption is 
critical: a model in which more than two KaiB and two 
KaiA dimers can bind also generates oscillations. 

4- The rate of spontaneous phosphorylation is lower than 
that of spontaneous dephosphorylation. 
The model includes spontaneous phosphorylation and de- 
phosphorylation of both active and inactive KaiC. Since 
KaiC reaches a low phosphorylation level in the absence 
of KaiA (and KaiB) 0, 0, [k|, the rate of spontaneous 
phosphorylation is lower than that of spontaneous de- 
phosphorylation. 
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FIG. 4: KaiC phosphorylation in the absence of KaiA and 
KaiB (C alone), in the presence of KaiA (C+A), in the pres- 
ence of KaiB (C+B) and in the presence of both KaiA and 
KaiB (C+A+B). For C+A, KaiC is initially fully unphospho- 
rylated; for C alone and C+B, KaiC is initially fully phospho- 
rylated (see also [i|). Parameters as in Fig. [5] 



FIG. 3: Sustained oscillations in the full model defined 
by Eqs. IrllllOl (A) The mean phosphorylation level p 
of KaiC shows a stable 24hr. rhythm. (B) Kai com- 
plexes. At t = 0, KaiC is fully unphosphorylated; 
[Co] = [C] T ; [A] = [A] T ; [B] = [B] T . KaiA then binds KaiC 
and stimulates phosphorylation. Next, the amount of KaiB- 
KaiC complex ([BC]) increases at high phosphorylation as 
KaiB binds to the inactive state of KaiC. Subsequently, KaiA 
is sequestered into a KaiA-KaiB-KaiC complex (ABC). The 
total concentrations equal those used in the in vitro experi- 
ments of 0: [C] T = 0.58/xM; [A] T = 1.75/uM; [B] T = 0.58AtM, 
corresponding to [A]t — [C]t and [B]t — 3[C]t- For other 
parameter values, see Table 1 of the Supporting Information. 



This model is described by the following reactions: 

.A ~ 
Cj <^ Ci 

bi 



Ci 
Ci 



A 



Kj 



AC, H f Ci+i + A 



2B <± B 2 Ci 



B 2 Ci + 2A z± A 2 B 2 Q 



C 



i+li 



^dpa 



B 2 Q ^ B 2 Ci + i, A 2 B 2 Ci 



A 2 B 2 C 



(6) 

(7) 
(8) 

(9) 

(10) 



As in Section [Til we assume the reaction rates are given 
by deterministic, mass-action kinetics. The most criti- 
cal parameters are the (de)phosphorylation rates. They 
have not been directly measured, but they are strongly 
constrained by the large number of quantitative in vitro 
experiments on the subsets of Kai proteins (see below). 
The model's predictions are much less sensitive to the 
remainder of its 39 parameters; for these we have simply 
chosen plausible values (see Supporting Information). 
KaiA + KaiB + KaiC Fig. [3K shows that our model 
produces sustained oscillations in KaiC phosphorylation 
when all three Kai proteins are present in the concentra- 



tions used in Q. Both the period and the amplitude of 
the oscillations agree well with those observed in [1, 0] . 
Fig. [3j3 shows the concentrations of complexes contain- 
ing KaiA and KaiC ([AC]), KaiB and KaiC ([BC]), and 
KaiA, KaiB, and KaiC ([ABC]), as a function of time. 
In the phosphorylation phase of the oscillations, KaiA 
binds to KaiC and stimulates its phosphorylation. At 
the top of the phosphorylation cycle, where KaiC hex- 
amers flip from the active to the inactive state, KaiA is 
released and KaiB binds to the inactive KaiC hexamers. 
The binding of KaiB stabilizes the inactive form of KaiC, 
preventing phosphorylation by KaiA. One critical role of 
KaiB is thus to allow the KaiC hexamers to enter the 
dephosphorylation phase of the cycle. 

Fig. [3)3 also shows that after [BC] has increased, [ABC] 
increases. This is because B 2 C4 — B 2 Ci can bind KaiA. 
This illustrates the second function of KaiB: KaiB that is 
bound to KaiC also sequesters KaiA. This leads to a form 
of the differential affinity mechanism at the end of the 
dephosphorylation phase of the cycle: The KaiC hexam- 
ers that are still in the inactive form — the laggards — will 
take away KaiA from those hexamers that have already 
flipped from the inactive to the active state — the front 
runners. This reduces the phosphorylation rate of the 
front runners, allowing the laggards to catch up. 

In our model, differential affinity acts at the bottom of 
the dephosphorylation phase of the cycle and throughout 
the phosphorylation phase. From the perspective of syn- 
chronizing the oscillations of the different hexamers, the 
ideal would be an ever-decreasing affinity between KaiA 
and KaiC, even as a given hexamer passes through the 
same sequence of states again and again. Thermodynam- 
ics, however, dictates that the affinity of KaiA for KaiC 
must increase somewhere in the cycle. In our model, this 
happens at the top of the inactive branch, where B 2 Cg 
and B 2 C5 do not bind KaiA, but B 2 C4 does have a high 
affinity for KaiA. 

KaiA + KaiC Fig. [4] shows that, in the presence of 
only KaiA, initially unphosphorylated KaiC reaches a 
phosphorylation level of about 90-95% after 6-8 hrs, in 
good quantitative agreement with experiment Q . In our 
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model, KaiC is biased towards the active state, and KaiA 
binding increases the stability of the active state even 
further. This explains the high steady-state phosphory- 
lation level when only KaiA is present. 
(KaiB +) KaiC Fig. 0] also shows that the phosphory- 
lation behavior of KaiC in the presence of KaiB is very 
similar to that of KaiC alone, as observed [H, [HI. Our 
model can explain this observation by assuming that the 
spontaneous dephosphorylation rate of the two KaiC con- 
formations is the same and is unaffected by KaiB bind- 
ing, which only stabilizes the inactive state with respect 
to the active one. 



B. Temperature compensation 

A striking feature of the in vitro oscillations of the Kai 
system is that they are temperature compensated [||. 
Specifically, as the temperature is increased from 25°C 
to 35°C, the period of the oscillations decreases by only 
10%. In general, the oscillation period of a network 
depends upon the rates of all the reactions in the sys- 
tem. In principle, one could try to achieve temper- 
ature compensation by balancing the temperature de- 
pendences of all of these rates [3fJ. We have adopted 
a different approach that is motivated by the fact that 
the (de) phosphorylation reactions are each individually 
temperature compensated 0] : The phosphorylation time 
courses of KaiC alone and of KaiC with KaiA change lit- 
tle between 25°C and 35°C. Indeed, the key idea of our 
approach is to construct the model so that the oscilla- 
tion period is determined by those rates that are known 
from experiment to be robust against temperature varia- 
tions while leaving it insensitive to the other rates, which 
might vary with temperature. 

A natural idea is to demand that the rates that 
can vary with temperature be much faster than the 
(de)phosphorylation rates, so that the period is domi- 
nated by the latter, which are temperature compensated. 
This leads to: 

5. All (un)binding rates and the flip rates fe and bo are 
much faster than the ( de )phosophorylation rates. 

Most conformational transitions are made at the top and 
bottom of the cycle; the period is thus less sensitive to 
flip rates other than fe and bo. 

Even when the (un)binding reactions between the 
Kai proteins are fast, however, the period can still de- 
pend upon the ratios of their rates — the dissociation 
constants — which will vary with temperature. The pe- 
riod becomes independent of the dissociation constants 
if all binding reactions go to completion. This occurs 
when the dissociation constants are much smaller than 
typical protein concentrations; in this limit, a change in 
the dissociation constants will have no appreciable effect 
on the fraction of bound proteins. We thus require: 

6. The affinities among the Kai proteins are high. 
KaiA, the least abundant of the three proteins, will then 
be almost entirely bound up in complexes with KaiB 
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FIG. 5: Temperature compensated oscillation period. The 
period of KaiC phosphorylation changes by 10% when the 
forward (f) or backward (pi) flip rates are changed by a factor 
25 (A) and by less than 5% when the dissociation constants 
for all KaiA and KaiB binding reactions are simultaneously 
changed by a factor 25 (B). Parameters as in Fig. 



and KaiC, in agreement with Q. As long as the relative 
magnitudes of the dissociation constants do not change 
with temperature, the composition of these complexes 
will moreover be unaffected. The phosphorylation rates, 
which depend on [AC,], are then robust to changes in 
temperature. Another important consequence of condi- 
tion 6 is that a proportional increase in all of the protein 
concentrations will have no effect on the oscillations, as 
has been observed @. 

Since no data on the temperature dependence of the 
dissociation constants and flip rates exists, we made the 
following estimate: We assumed that both the binding 
energies and the energy barriers for the conformational 
transitions are at most 50fceT. If the temperature is 
changed from 25° to 35°, the dissociation constants and 
flip rates can then change by about an order of magni- 
tude. To test whether our model is robust against such 
perturbations, we have varied both dissociation constants 
and flip rates by a factor five in each direction. Fig. [5] 
shows that our model withstands these trials: The pe- 
riod varies by about 5 — 10%, in very good agreement 
with experiment This is strong evidence that con- 
ditions 5 and 6, together with temperature-compensated 
(de)phosphorylation rates, are sufficient for temperature- 
compensated oscillations. 



KaiC dynamics as a function of the KaiA and 
KaiB concentration 



Fig.[S]shows the behavior of our model as a function of 
the total KaiA and KaiB concentrations [A]t and [B]t- 
For [A]t < 0.5[C]t, the system exhibits no oscillations. 
At around [A]t = 0.5[C]t, the system starts to oscil- 
late via a supercritical Hopf bifurcation with a period 
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FIG. 6: KaiC oscillations as a function of KaiA and KaiB 
concentration. (A) Period and (B) amplitude of oscillations 
in KaiC phosphorylation as a function of the concentration of 
KaiA and KaiB. The dashed curve shows the location of the 
supercritical Hopf bifurcation that gives birth to the oscilla- 
tions, and the color scales give period in hours and amplitude 
of p oscillation. Note the small region of bistability (solid line; 
see also Supporting Information) at low [A]t and [B]t- The 
remaining parameters are as in Fig. [3] (C) KaiC oscillations 
as a function of KaiA concentration. Results are shown for 
[B] T = 3[C] T and [A] T = 0.2[C] T (green), 0.6[C] T (blue), [C] T 
(black) and 1.4[C]t (red). (D) KaiC oscillations as a function 
of KaiB concentration. Results are shown for [A]t = [C]t and 
[B] T = 1.2[C] T (purple), 2.1[C] T (yellow) and 3[C] T (black). 



of about 35 hrs (see Supporting Information for details 
on the bifurcation analysis). As the KaiA concentration 
is increased, the period monotonically decreases. In con- 
trast, the amplitude first increases to reach a maximum 
at around [A]t = 0.85[C]t, then decreases until oscilla- 
tions disappear at around [A]t = 1.25[C]t- The dynam- 
ics as a function of the KaiB concentration are markedly 
different. Fig. [6] shows that a minimum KaiB concen- 
tration of about [B]t = [C]t is needed to sustain oscilla- 
tions. Above that threshold, neither the period nor the 
amplitude depend strongly on [B]t- 

The different effects of varying [A]t and [B]t can be 
understood from the different roles the two dimers play 
in our model. KaiA stimulates the phosphorylation of 
KaiC. If the total KaiA concentration is very low, the 
phosphorylation rate will thus be so slow that it is coun- 
terbalanced by the spontaneous dephosphorylation rate. 
If, on the other hand, the total KaiA concentration is 
very high, the mechanism of differential affinity no longer 
functions, because it relies on competition for a limited 
amount of KaiA. The function of KaiB is to stabilize in- 
active KaiC and to sequester KaiA. As long as enough 
KaiB is available to perform these functions, the period 
and amplitude will not depend upon the KaiB concen- 
tration. 

Interestingly, the very recent experiments of Q give 



strong support for our model. In particular, these exper- 
iments show that when the KaiA and KaiB concentra- 
tions are reduced from their standard values by a factor 
of four and three, respectively, all oscillations cease in 
very good agreement with our results. We further pre- 
dict that there is an upper bound on the KaiA concentra- 
tion, but not on the KaiB concentration, for oscillations 
to exist. Moreover, while the amplitude and the period 
of the oscillations do not depend in our model on [B]t, 
they do depend in a very specific manner on [A]x- These 
dependences on the KaiA and KaiB concentrations are 
direct consequences of the basic roles of these proteins in 
our model. They thus represent some of our most robust 
and important predictions. 



IV. DISCUSSION 

We have presented an allosteric model of KaiC phos- 
phorylation that can describe a wealth of experimental 
data on the Kai system. Its foundation is the assumption 
that each KaiC hcxamcr can exist in two distinct confor- 
mational states, an active one in which it tends to be 
phosphorylatcd and an inactive one in which it tends to 
be dephosphorylatcd. Because of the interplay between 
nucleotide binding, which favors the active state, and 
phosphorylation, which favors the inactive state, each 
individual hexamer will repetitively gain and lose phos- 
phate groups. However, if macroscopic oscillations are to 
be observed, the phosphorylation cycles of the individual 
hexamers must be synchronized. We introduced a novel 
mechanism, called differential affinity, which, in contrast 
to some previous models [2(J [H|, allows for synchroni- 
sation even in the absence of direct interactions between 
hexamers. The key idea is that while all KaiC hexam- 
ers compete to bind KaiA, which stimulates phosphory- 
lation, the laggards in the cycle are continuously being 
favored in the competition. This mechanism is most ef- 
fective when KaiB and KaiC bind KaiA very strongly. 
It is also precisely in this limit that the oscillation pe- 
riod becomes insensitive to changes in the Kai proteins' 
affinities for each other. Differential affinity and temper- 
ature compensation are thus intimately connected. The 
mechanism of driving two-body reactions to saturation 
is, however, more general; it could, for instance, be used 
to make temporal programs of g ene expression robust 
against temperature variation [ldj | . 

In S. elongatus, the concentration of KaiA dimers is 
less than 10% of that of KaiC hexamers Our model 
predicts that in this regime the in vitro oscillations of [J 
disappear. The very recent experiments of Q support 
this prediction: they unambiguously demonstrate that 
in vitro the oscillations cease to exist if the concentra- 
tion of the KaiA dimers is less than 25% of that of the 
KaiC hexamers. Clearly, in vivo other processes are at 
work. It is known, for instance, that both the subcellular 
localization of the Kai proteins [l3[ and KaiC's role as 
a transcriptional repressor Q affect circadian rhythms, 
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as do other clock proteins such as SasA It is tempt- 
ing to speculate, however, that these additional effects 
merely shift the phase boundaries of the model presented 
here without changing its basic mechanism. One could 
imagine, for example, that a combination of KaiB lo- 
calization to the cell membrane and competitive binding 
by molecules like SasA could reduce the number of sites 
available to sequester KaiA, thus allowing the oscillator 
to function at lower KaiA concentrations. 

Finally, our model makes a number of predictions that 
could be verified experimentally. One clear prediction is 
that KaiC can exist in two distinct conformational states. 
Moreover, our model suggests that KaiC binds KaiA and 



KaiB very strongly, with dissociation constants that de- 
pend upon the conformational state and phosphorylation 
level of the KaiC hexamer. But perhaps the strongest 
test of our model concerns the KaiC oscillation dynam- 
ics as a function of the KaiA and KaiB concentrations 
(see Fig. [6]): We predict that the oscillations will disap- 
pear when the KaiA concentration is increased, but not 
when the KaiB concentration is increased. 

We thank Daan Frcnkcl and Martin Howard for a crit- 
ical reading of the manuscript and the Aspen Center for 
Physics for its hospitality to D.K.L. early in this project. 
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